BABAR-PUB-03/021 
SLAC-PUB-10132 



Evidence for the Rare Decay B — > K*£ + £ and Measurement of the B — > K£ + £ 
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We present evidence for the flavor-changing neutral current decay B —> K*£ + £~ and a mea- 
surement of the branching fraction for the related process B — > K£ + £~, where £ + £~ is either an 
e + e~ or pair. These decays are highly suppressed in the Standard Model, and they are 
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sensitive to contributions from new particles in the intermediate state. The data sample com- 
prises 123 x 10 6 T(4S) -> BB decays collected with the BABAR detector at the PEP-II e + e" 
storage ring. Averaging over isospin and lepton flavor, we obtain the branching fractions 

B{B -> Ki+r) = (0.65±g;if ±0.04) x KT 6 and B(B -> K*l + i~) = (0.88±g:|| ±0.10) x KT 6 , where 
the uncertainties are statistical and systematic, respectively. The significance of the B — > K£ + £~ 
signal is over 8<r, while for B — * K*£ + £~ it is 3.3<r. 

PACS numbers: 13.25.Hw, 13.20.He 



Rare decays of B mesons that involve loop diagrams 
in the Standard Model (SM) provide a promising means 
to search for effects beyond the SM yj. The decays 
B -> K£+£~ and B -> K*£ + £~, where £± are charged 
lcptons and K* is the if* (892) meson, result from one- 
loop processes that transform the 6-quark in the initial- 
state B meson into an s-quark in the final-state K^*' 
meson. 

In the SM, three amplitudes contribute at leading or- 
der: an electromagnetic (EM) penguin, a Z penguin, and 
a W + W~ box diagram. The penguin diagrams involve 
the emission and absorption of a W boson. The pres- 
ence of three SM electroweak amplitudes makes B — > 
K (*>£ + £~ more complex than B — > if*7, which proceeds 
solely through an EM penguin. 

Because of their loop structure, these decays are highly 
suppressed, with SM branching fractions expected to be 
roughly 0.5 x 10~ 6 for B — » K£ + £~ and about three times 
that for the B -> K*i+(r modes E 0. Due to the 
complexity of strong interaction effects, however, theoret- 
ical uncertainties on the rates are currently at least 35% 
(Ali et al. 1]). Both B -» K*e+p- and B -> K*n + n~ 
receive a contribution from the pole in the EM penguin 
amplitude at q 2 = rn 2 +g _ — 0, but the enhancement in 
the electron mode is larger. An important consequence of 
the loop structure of these decays is that their branching 
fractions and kinematic distributions can be significantly 
affected by the presence of new particles, such as those 
predicted in models based on supersymmetry 0. 

Recently, substantial progress has been made in ex- 
perimental studies of these decays. The Belle collabora- 
tion has observed B — > K£ + £~, as well as the inclusive 
B — > X s £ + £~ decay 0. Limits on these and similar 
modes have been set by BABAR 0, Belle 0, CLEO 0, 
and CDF 0. Our new measurements are based on a 
data sample six times larger than that used for our pre- 
viously published results. We study eight final states: 
B+ -> K+£+£r, B° -> K%£+£-, B+ -> K*+£+£~, and 
B° -> K*°£ + £~, where K*° -> K + Tr~, K*+ -> K s ir + , 
Kg — > tt + it~ , and £ is either an e or fx. Throughout this 
paper, charge-conjugate modes are implied. 

We analyze data collected with the BABAR detec- 
tor at the PEP-II storage ring at the Stanford Lin- 
car Accelerator Center. The data sample comprises 
113.1 fb _1 recorded on the T(4S) resonance, yielding 
(122.9 ± 1.4) x 10 6 BB decays, and an off-resonance sam- 
ple of 12.0 fb _1 used to study continuum background. 



We select events that include two oppositely charged 
lcptons (e + e~, /.t + /i~), a kaon (either or Kg), and, 
for the B — > K*£ + t~ modes, a tt^ 1 that combines with 
a kaon to form a K* candidate. Electrons are identi- 
fied primarily in the CsI(Tl) electromagnetic calorimeter, 
while muons are identified by their penetration through 
iron plates of the magnet flux return. Electron (muon) 
candidates are required to satisfy p > 0.5 (1.0) GeV/c. 
Bremsstrahlung photons from electrons are recovered by 
combining an electron candidate with up to one photon 
with Ej > 30 MeV in a small angular region around 
the initial electron direction. Photon conversions and ir° 
Dalitz decays are removed by vetoing all low-mass e + e~ 
pairs, except in B — > K*e + e~ modes, where we preserve 
acceptance at low mass by retaining pairs that intersect 
inside the beam pipe. 

K ± candidates are tracks with dE/dx and Cherenkov 
angle consistent with a kaon. tt^ 1 candidates are tracks 
that do not satisfy the if selection. K s candidates are 
reconstructed from two oppositely charged tracks with 
an invariant mass consistent with the K s mass and a 
common vertex displaced from the primary vertex by at 
least 1 mm. 

True B signal decays produce narrow peaks in the dis- 
tributions of two kinematic variables, which can be fit- 
ted to extract the signal and background yields. For a 
candidate system of B daughter particles with masses 
m% and three-momenta p* in the T(AS) center-of-mass 
(CM) frame, we define m E s = ^/i?,* 2 — | J2i P*l 2 anc ^ 
AE = J2i \/ m i + Pi — -®b! where E£ is the beam energy 
in the CM frame. For signal events, the mES distribu- 
tion peaks at the B meson mass with resolution a w 
2.5 MeV/c 2 , and the AE distribution peaks near zero, 
with a typical width a m 20 MeV. In B -> K£+£~ chan- 
nels, we perform a two-dimensional unbinned maximum- 
likelihood fit to the distribution of toes and AE in the 
region to es > 5.2 GeV/c 2 and \AE\ < 0.25 GeV. In 
B — > K*£ + £~ decays, we perform a three-dimensional fit 
to TOesi and the kaon-pion invariant mass in the 

region 0.7 < m K7I < 1.1 GeV/c 2 . 

Backgrounds arise from three main sources: random 
combinations of particles from qq events produced in 
the continuum, random combinations of particles from 
T(4S) — * BB decays, and B decays to topologies simi- 
lar to the signal modes. The first two ("combinatorial") 
backgrounds typically arise from pairs of semileptonic de- 
cays and produce broad distributions in toes arid AE 
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compared to the signal. The third source arises from 
modes such as B — > J/ipK^-*' (with J /if) — > ^ + £~) or 
£? — > K^ttit (with pions misidentified as muons), which 
have shapes similar to the signal. All selection criteria 
are optimized with GEANT4 9] simulated data or with 
data samples outside the full fit region. 

We suppress combinatorial background from contin- 
uum processes using a Fisher discriminant which 
is a linear combination of variables with coefficients op- 
timized to distinguish between signal and background. 
The variables (defined in the CM frame) are (1) the ratio 
of second- to zeroth-order Fox- Wolfram moments for 
the event, computed using all charged tracks and neutral 
energy clusters; (2) the angle between the thrust axis of 
the B candidate and that of the remaining particles in 
the event; (3) the production angle 8 b of the B candidate 
with respect to the beam axis; and (4) the masses of Kl 
pairs with charge correlation consistent with D decay. 

We suppress combinatorial backgrounds from BB 
events using a likelihood function constructed from (1) 
the missing energy of the event, computed from all 
charged tracks and neutral energy clusters; (2) the vertex 
fit probability of all tracks from the B candidate; (3) the 
vertex fit probability of the two leptons; and (4) the an- 
gle 9b- Missing energy provides the strongest suppression 
of combinatorial BB background events, which typically 
contain neutrinos from two scmilcptonic decays. 

The most prominent backgrounds that peak in niES 
and AE are B decays to charmonium: B — > J/tpK^*' 
(with J/ip — > l + £~) and analogous B decays to i/j(2S). 
We exclude dilepton pairs consistent with the J/tp 
(2.90 < m e+e - < 3.20 GeV/c 2 and 3.00 < m M+li - < 
3.20 GeV/c 2 ) or with the ^(25) (3.60 < m l+t - < 
3.75 GeV/c 2 ). This veto is also applied to m e+e - com- 
puted without bremsstrahlung photon recovery. When 
a lepton radiates or is mismeasured, mi+i- can shift 
away from the charmonium mass, while AE shifts in 
a correlated manner. The veto region is extended in 
the {m,£+i-,AE) plane to account for this correlation, 
removing nearly all charmonium events and simplifying 
the description of the background in the fit. Because 
the charmonium events removed by these vetoes are so 
similar to signal events, these modes provide extensive 
control samples (about 5200 events in all) for studying 
signal shapes, selection efficiencies, and systematic er- 
rors. Outside the charmonium veto regions, the signal 
efficiency is similar over the full q 2 range of each mode. 

In muon modes, where the probability for a hadron 
to be misidentified as a muon can be as high as a few 
percent, background from the decay B~ — > D°ir~ with 
D° -> R-tt+_ or D° -> K*-ir+, or from B° -> D+ir- 
with D + — ► K*°tt + , is significant. These events are sup- 
pressed by vetoing events where the K^fi kinematics are 
consistent with those of a hadronic D decay. 

We estimate the residual peaking background from 
measurements in the data, supplemented in some cases 



TABLE I: Results from the fits to B -» K^l+T modes. 
The columns are, from left to right: fitted signal yield; the 
signal efficiency, e (not including the branching fractions for 
K* and K° decays); the fractional systematic error on the se- 
lection efficiency, AB € /B e ; the systematic error from the fit, 
ABfi t ; and the branching fraction central value (B) with its 
statistical and total systematic uncertainties. For the branch- 
ing fractions averaged over different channels (lower part of 
table), simultaneous, constrained fits are performed to ex- 
tract an efficiency-corrected signal yield that averages over 
the included channels. The modes with significance > 3a 
are K + e+e~ (8.4a), (4. Iff), Ke+e~ (7.8a), K£+£~ 

{8 Act), and K* 1+ f (3.3a). 





Signal 


e AB C /B 


A£ flt 


B 


Mode 


yield 


(96) (96) 


(to- 6 ) 


(10-6) 



K+e+e- 

K+n+ti- 

K°e+e- 

K /i+fi~ 

K*°e+e- 

K*°ii+fi' 

K'+e+e- 



24.7 
0.7 

-1 
5.9 

12.4 
4.5 
0.6 
4.2 



+ 5.9 
5.2 
2.0 
1.2 
+2.0 
-1.4 
+3.0 
2.3 
+6.3 
5.2 
+4.1 
3.0 
+3.8 
2.5 
+3.5 
2.4 



19.2 

8.5 
20.1 

8.6 
13.6 

6.4 
10.2 

4.8 



±6.3 
±7.6 
±8.4 
±8.8 
±7.6 
±10.1 
±10.7 
±12.7 



±0.02 
±0.02 
±0.08 
±0.02 
±0.08 
±0.07 
±0.28 
±0.15 



1.05±° 

0.07±o 
i+o 

o 

+ 

+0 

+ 


+1 



+2 
1 



-0.21 
1.63 
1.11 
0.86 
0.20 
3.07 



.22 ±0 - 07 
;^±0.02 

■^±0.08 

;g±o.i4 

if ±o.n 

^±0.28 

'.78 ±0 - 42 



Ke+e- 

Kt+tr 

K*e+e~ 
K*l+l~ 



91 
55 
80 
121 
156 
108 



+22 
1!) 
+29 
23 
17 
-15 
+ 01 
51 
+ 94 
75 
+41 
36 



±6.5 
±7.4 
±6.4 
±7.8 
±10.1 
±8.1 



±0.02 
±0.01 
±0.01 
±0.08 
±0.09 
±0.07 



0.74l°;^±0.05 
s+o.23 ±ao4 



0.45Z 
0.65j 
0.98 11 



-0.19- 

5l°i^±0.04 
3±°i°±0.11 
1.27t^6±0.16 

d+0.'33_ 



0.88] 



0.2!) 



±0.10 



by simulation studies. Events from B — > K^rnr, B — > 
K^Ktt, and B — > K^KK are highly suppressed by the 
particle identification criteria. These backgrounds are es- 
timated from control samples to be 0.19 ± 0.11 events 
per channel averaged over muon modes and less than 
0.01 events per channel in electron modes. After the ve- 
toes on B -> J/ipKW and B -> tp(2S)K^ decays, the 
remaining peaking background is estimated from simula- 
tion to be 0.17 ± 0.07 events per channel averaged over 
B -> K*£+£- modes, and it is negligible in B -> K£+£~ 
modes. The background from B — > K*j (with photon 
conversion in the detector) is determined from simulation 
to be 0.48±0.16 events in B° -> K*°e+er and 0.09±0.04 
events in B + — * K* + e + e~. 

The signal shapes are parameterized with a Gaussian 
core for ttt-es and a double Gaussian core for AE. Both 
the mES ^md AE shapes include a radiative tail, which 
accounts for the effects of bremsstrahlung. The ibes 
shape parameters are assumed to have AE dependence 
Co + C2{AE) 2 . All signal shape parameters are fixed from 
signal simulation, except for the mean and width param- 
eters in mES (co only) and AE, which are fixed to values 
from charmonium data control samples. 
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The background is modeled as the sum of three terms: 
(1) a combinatorial background shape with floating nor- 
malization, written as the product of an ARGUS func- 
tion [l2T ] in wes, an exponential in AE, and the prod- 
uct of ^/rriKiT — ttik — ttItt and a quadratic function of 
mK-n for the K* modes; (2) a peaking background con- 
tribution, with the same shape as the signal, but with 
normalization fixed to measured peaking backgrounds; 
and (3) terms with floating normalization to describe 
(a) background in B -> K£+£~ (B -> K*£+£') from 
B -> K*e + fr (B ->■ K*ir£ + £-) events with a lost pion, 
and (b) background in B -> K*£+£- from B -> K£+£~ 
events with a randomly added pion. In the K* modes, 
we allow an additional background (4) that uses our com- 
binatorial shape in toes and AE, but peaks in rrijcir at 
the K* mass. Because the normalizations for terms (1), 
(3), and (4) are floating, as are the combinatorial back- 
ground shape parameters, much of the uncertainty in the 
background is propagated into the statistical uncertainty 
on the signal yield obtained from the fit. 

Table [I] lists signal yields and branching fractions for 
each mode. The relative systematic uncertainties on the 
efficiency, AB e /B, arise from charged-particle tracking 
(1.0% per lepton, 1.7% per charged hadron), particle 
identification (1.1% per electron, 1.6% per muon, 0.9% 
per pion, 0.9% per kaon), the continuum suppression cut 
[(0.8-2.8)%], the BB suppression cut [(1.4-5.0)%], K° s 
selection (3.8%), signal simulation statistics [(0.7-1.4)%], 
theoretical model dependence of the efficiency [(4-7)%, 
depending on the mode], and the number of BB events 
(1.1%). Uncertainties on efficiencies due to model de- 
pendence of form factors are taken to be the full range 
of variation from a set of models . 

The systematic uncertainties on the fit yields, AB^t, 
arise from three sources: uncertainties in the parameters 
describing the signal shapes, possible correlation between 
ibes and AE in the combinatorial background shape, 
and uncertainties in the peaking backgrounds. The un- 
certainties in the means and widths of the signal shapes 
are obtained by comparing data and simulation for the 
charmonium control samples. For modes with electrons, 
we also vary the fraction of signal events in the tail of 
the AE distribution. To evaluate sensitivity to the back- 
ground parameterization, we allow additional parameters 
and a correlation between toes and AE. 

Table[I]also lists results from simultaneous fits to com- 
binations of B — > K£ + £~ modes and combinations of 
B — > K*£ + £~ modes, where the relative branching frac- 
tions for the contributing modes are constrained. B° 
and B + production rates are constrained to be equal, 
and the ratio of their total widths is constrained to be 
1.085 ± 0.017 0]. All branching fractions from simulta- 
neous fits are expressed in terms of the B° total width. 
The projections of the fit on toes and AE are shown in 
Fig. n for the simultaneous fit to the B — > K£ + £~ chan- 
nels. We assume that all four B — > K£ + £~ modes have 




5.2 5.22 5.24 5.26 5.28 -0.2 -0.1 0.1 0.2 

m ES (GeV/c 2 ) A E (GeV) 

FIG. 1: Distributions of the fit variables in K£ + £~ data 
(points), compared with projections of the simultaneous fit 
(curves): (a) ttibs distribution after requiring —0.11 < AE < 
0.05 GeV and (b) AE distribution after requiring |mES — 
ms| < 6.6MeV/c 2 (2.6a). The solid curve is the sum of all fit 
components, including signal; the dashed curve is the sum of 
all background components. 

equal partial widths. A signal is evident at the B mass 
in toes and at AE = 0. Figure shows projections of 
the simultaneous fit to all B — ► K*£ + £~ modes. Here, 
the partial width ratio of electron and muon modes is 
constrained to be T(B -> K*e + e^)/T(B -> K*fi + fi~) = 
1.33 from the model of Ali el al. [l|. Our simultaneous 
fit result is expressed as a B° — > K*°[i + ij,~ branching 
fraction. 

The significance of the B — > K£ + £~ signal from the 
simultaneous fit is ~ 8<r, computed as \/2A\ogC, where 
Alog£ is the likelihood difference between the best fit 
and the null-signal hypothesis. We account for system- 
atic uncertainties in the significance by simultaneously 
including all effects that individually lower the fit yields 
prior to computing the change in likelihood. The signifi- 
cance of the B — > K*£ + £~ signal, including all systematic 
uncertainties, is 3.3cr (3.8cr not including them). 

In summary, we have observed signals for B — > K£ + £~ , 
averaged over lepton type (e + e _ and and B 

charge, and we have obtained the first evidence for 
B — > K*£ + £~ , similarly averaged. We obtain 

B(B -> K£+£-) = (0.651°;^ ± 0.04) x 10~ 6 , 
B(B -► K*l+r) = (0.88t£! ± 0.10) x 10~ 6 , 

where the first error is statistical and the second is sys- 
tematic. Our branching fraction for B — > K£ + £~ is 
slightly higher than our previous limit 0.51 x 10 -6 (90% 
confidence level) Q and is in agreement with the Belle 
result (0.75±g;|f ± 0.09) x lO" 6 4]. Our B -► K*£+£~ 
branching fraction is consistent with previous 90% confi- 
dence level limits from BABAR (< 3.1 x 10~ 6 for K*t + l~) 
and Belle (< 3.1 x 10~ 6 for K*^ + ^-) 4]. These re- 
sults are consistent with the range of predictions based 
on the Standard Model [HSU. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
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FIG. 2: Distributions of the fit variables in K*£ + £~ data 
(points), compared with projections of the simultaneous fit 
(curves): (a) jties after requiring —0.11 < AE < 0.05 GeV 
and 0.817 < m Kl , < 0.967 GeV/c 2 , (b) AE after requir- 
ing |m E s - m B \ < 6.6MeV/c 2 {2.6a), 0.817 < mic* < 
0.967 GeV/c 2 , and (c) mx» after requiring |tties — tub\ < 
6.6MeV/c 2 and -0.11 < AE < 0.05 GeV. The solid curve 
is the sum of all fit components, including signal; the dashed 
curve is the sum of all background components. 
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